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STRUCTURE IN THE ROTATION MEASURE SKY 
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ABSTRACT 

An analys i s of s tructure in rotation measure (RM) across the sky based on the RM catalog of 
I Tavlor et al.l (|2009t ) is presented. Several resolved RM structures are identified with structure in the 
local ISM, including radio loops I, II, and III, the Gum nebula, and the Orion-Eridanus super bubble. 
Structure functions (SFs) of RM are presented for selected areas, and maps of SF amplitude and 
slope across the sky are compared with Ha intensity and diffuse polarized intensity. RM variance 
on an angular scale of 1° is correlated with length of the line of sight through the Galaxy, with a 
contribution from local structures. The slope of the SFs is less concentrated to the Galactic plane 
and less correlated with length of the line of sight through the Galaxy, suggesting a more local origin 
for RM structure on angular scales ~ 10°. The RM variance is a factor ~ 2 higher towards the SGP 
than towards the NGP, reflecting a more wide-spread asymmetry between the northern and southern 
Galactic hemispheres. Depolarization of diffuse Galactic synchrotron emission at latitudes < 30° can 
be explained largely by Faraday dispersion related to small-scale variance in RM, but the errors allow 
a significant contribution from differential Faraday rotation along the line of sight. 
Subject headings: ISM: magnetic fields, ISM: structure, radio continuum: ISM 



1. INTRODUCTION 

The magnetic field of the Milky Way Galaxy is the 
only galactic magnetic field that can be studied over ~ 
5 orders of magnitude, on length scales less than a parc- 
scc to several kiloparsec. The structure of this magnetic 
field is revealed by different observing techniques that are 
usually sensitive to either the line-of-sight component of 
the magnetic field (Faraday rotation, Zeeman splitting), 
or the component of the magnetic field perpendicular to 
the line of sight (stellar polarization, polarization of radio 
synchrotron emission). Zeeman splitting provides a di- 
rect in-situ measurement of the strength of the magnetic 
field, while Faraday rotation and stellar polarization are 
propagation effects that require additional information 
on the electron density or the extinction column of dust 
along the line of sight. Galactic synchrotron emission 
at frequencies of a few GHz suffers from strong depolar- 
ization from Faraday rotation in the Galactic disk (e.g. 
Sokoloff ct al. I fliMfGaensler et al.ll200ll : lUvaniker et al.l 
20031). Even at high, frequencies (e.g. 20 GHz WMAP; 
Page et all I2007T ) where Faraday rotation is not signifi- 
cant, the summation of polarized emission with different 
polarization angles over a long line of sight can be a sig- 
nificant source of depolarization. 

The large-scale structure of the Galactic magnetic field 
has been revealed mostly by Faraday rotation of polar- 
ized emissio n from pulsars and distant active galac t ic nu- 
clei (AGN) (iSimard-Normandin fe Kronberd|l979l Il980l; 
Han et al.l 12006b iBrown et al.l I2001L I2003L 120071: iValled 
2008HMen etaLll2008t iTavlor et al.ll2009l among others). 
The angle change A9 from Faraday rotation of polarized 
emission of a compact radio source with negligible inter- 
nal Faraday rotation is given by 



A9 = A 2 0.8I2 



observer 



density, B is the magnetic field vector, and dl is a line of 
sight segment pointing from the source to the observer 
so that RM has a positive sign if there is a net magnetic 
field component along the line of sight in the direction of 
the observer. The angle change A9 can be determined 
from observations of the polarization angle of a radio 
source as a function of wavelength. 

The electron density along the line of sight must be 
determined independently to derive the strength of the 
magnetic field. For the large-scale magnetic field of the 
Galaxy, a model of the distribution of the electron den- 
sity n e along the line of si ght has been derived f rom pul- 
sar dispersion measures (ICordes fc Laziol[200l . How- 
ever, for smaller length scales, an electron density distri- 
bution must be proposed as an assumption that is largely 
unconstrained by current observations, as significant an- 
gle changes can occur in an ionize d region that cannot b e 
detected by its thermal emission ()Uvaniker et aLl feOOS). 
Correlation between fluctuations in electron density and 
magnetic field strength can bias the RM upward or down- 
ward, depending whether there is a positive correlation 
or an anti-correlation respective ly between therm al elec- 
tron density and field strength (Be ck et al . 20031) . 

Structure in the Galactic magnetic field on smaller 
scales has been apparent in polarization data of stars 
within 1 kpc fe.g.lAxon fc Ellis 1976; Mat hewson fc Ford 
19701: ILerovlll999T IHeilesI I200ct IBerdvugin fc Teerikorpi 
2002t IBerdvugin et al.l 12004 and references therein), 
and p olarization of diffuse Galac t ic synchrotron emission 
(e.g. iDuncan et al I |1997l 119991: lUvaniker et al"~l 119991: 
IWolleben et al.l 120061: iPage et al 
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.B • dl = X 2 RM (1) 



where A if the observing wavelength, n e is the electron 



Gacnsler et al. 

2007t ISchnitzeler et al.ll2007l: IWolleben ll2007tlGao et al 
120101 Landecker et al., in prep.). Diffuse polarized 
emission has the advantage that it provides contin- 
uous sampling of a foreground Faraday screen, but 
the signal is usually a complicated integral of syn- 
chrotron emissivity and Fa raday rotation over the vol- 
ume probed by the beam (Sokoloff ct al. 1998). Equa- 
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tion [T] does not apply in this case, but imaging po- 
larized intensity and polarization angle, and RM syn- 
thesis (|Burnl 119661: iBrentiens fc De Bruvnl 120051) have 
revealed structur e s in diffuse polarize d emission (e.g. 
Gray et all 119981: IGaensler et al.1 120011: lUvaniker et al.1 
2003HSchnitzeler et al.ll2007t iDe Bruvn et al.ll2009f). 

RM structure functions (SFs) dSimonetti. et al. 119841: 
iSimonetti fc Cordeslfl98llMinter fc Spanglerlll996!) have 
been used to obtain a statistical description of RM vari- 
ance as a function of angular scale. At high Galactic 
latitudes the varianc e in RM was found to be indepen- 
dent of angular scale (jSimonetti. et al.|[l984t ISun fc Hani 
I2004D . At lower latitudes, RM varia nce is often observed 
to increase wit h ang u lar scale (ISimonetti fc Cordesl 
119861: iSimonettil 119921 ISun fc Hani 120041) . sometimes 
with a break around an angular scale of 1° above 
which the variance remains constant or increases at a 
slower rate (lHaverkorn et al.ll2003l 120061 : lHaverkornl2007t 
lHaverkorn et al.ll2008HRov et al.ll2008D . The break in the 
RM SF is usually interpreted as an outer scale in the tur- 
bulence that gives rise to the RM variance. The smallest 
scale that can be probed depends on the density of RM 
measurements, and angular resolution of the radio data 
(ISimonetti fc Cordesl 1 19861: lLeahvlll987D . IGaensler et al.l 
( 2005) and lHaverkorn et al.l (|2008D used depolarization 
of point sources as a measure of RM variations on the 
angular scale of extragalactic radio sources. 

A database of 37 ,543 RMs that covers 80% of the sky 
(jTavlor et al.1 [2009) was construct ed recently from th e 
NRAO VLA Sky Survey (NVSS; ICondon et al~lll998l) . 
using the two frequency bands of the NVSS separated by 
70 MHz (center to center). The NVSS RMs are reliable 
to \RM\ < 520 rad m~ 2 . Higher RMs are subject to 
ambiguities because it was not possible to derive a unique 
RM using the available constraints and noise in the data. 
This occurs mainly at low Galactic latitude (|6| < 5°). 

2. LARGE-SCALE RM STRUCTURES AND THE LOCAL ISM 

iTavlor et al.l (|2009f) showed the signature of the large- 
scale magnetic field in the solar neighborhood in these 
data including a reversal of the sign of RM across the 
Galactic plane consistent with a quadrupole magnetic 
field geometry. Assuming a symmetric electron density, 
they found that the magnetic field strength south of the 
Galactic plane is approximately two times the magnetic 
field strength in th e north. Towards the Galactic poles, 
ITavlor et all (|2009) found a positive RM in the north 
and in the south, also stronger in the south. ? found 
a significant mean RM in the region around the South 
Galactic Pole (GPP), but not around the North Galactic 
Pole (NGP). 

Here we discuss RM structures on angular scales < 90° 
and their relation the local interstellar medium. Fig- 
ures [U [21 and [3] show the R M database as cir cles on a 
background of Ha intensity dFinkbeinerl 12003ft , diffus e 
polarized emission at 1.4 G Hz dWolleben et al l [20061) . 
and diffuse X-ray emission (jSnowden et al.l 119981) . re- 
spectively. Figure [4] shows a different representation of 
the RM sky towards the anticenter, excluding data with 
\RM\ < 25 rad m~ 2 , and showing all circles at the same 
size. This representation emphasizes RM structure with 
amplitudes of a few tens rad m -2 , but also boundaries 
where the sign of RM changes. 

Three large regions previously identified by 



Simard-Normandin fc Kronberg (1980) as regions 
A, B, and C are responsible for most of the structure 
on angular scales > 30°. These structures are also 
represe nted in other all-sky RM images with resolution 
~ 20° (iFrick et al. I l200lHJohnston^Hollitt et all 120041: 
iDineen fc Coles ll2005HXu et al.ll2006l ). The denser sam- 
pling of the present data reveals structure in more detail, 
and we adop t somewhat different boundaries for regi ons 
A and C of iSimard-Normandin fc K ronberg ()1980D as 
discussed below. The boundaries adopted here are listed 
in Table [IJ and we will refer to region A' and region C 
instead. Region B in Simard-Normandin & Kronberg 
(1980) refers to the Gum nebula, and will be discussed 
as such. 

Region A is a large area of negati ve RM for which 
ISimard-Normandin fc Kronberg! (jl980f ) listed the bound- 
aries as 60° < / < 140°, -40° < b < 10°. With 
the superior sampling of the new data, we see sub- 
stantial substructure in this region. The most nega- 
tive RMs occur at intermediate latitudes in the region 
80° < I < 150°, -40° < b < -20°. Distinctly smaller 
RM amplitudes are found in the area 100° < I < 150°, 
—20° < b < — 5°, except for an area with higher RM am- 
plitudes in the longitude range 120° < I < 130° (See also 
Figured]). The mean RM in the area 100° < I < 120°, 
-40° < b < -20° is -67.7 rad m~ 2 , with a standard de- 
viation of 38.5 rad m~ 2 . The mean Ha intensity in this 
area is 1.8 R. The mean RM in the area 100° < I < 120°, 
-20° < b < -5° is -32.8 rad m~ 2 with a standard 
deviation of 31.9 rad m -2 . The mean background Ha 
intensity in this area, avoiding bright H n regions is 
5 R. The section of region A with higher RM ampli- 
tude is therefore separated from the Galactic plane by 
a region with RM amplitudes a factor ~ 2 lower than 
those at more negative latitudes, even though the Ha in- 
tensity is a factor ~ 2.5 higher than at m o re ne gative 
latitudes. ISimard-Normandin fc Kronberg! ([19801 ) sug- 
gested that region A extends across the Galactic plane, 
to b = +10°, but this appears to be based on only two 
sources in their Figure 1, and there is no indication of 
an exte nsion to positive latitudes in thei r Figu re 2. The 
data of ISimard-Normandin fc Kronberg] (]1980l ) sampled 
the intermediate region with low RM amplitudes very 
sparsely. 

Region A is bordered on the high-longitude side by 
a line where RM nearly vanis hes, and then reverses 
sign to become positive (e.g. IXu et al.l 12006!) . The 
boundary with near-zero RM runs approximately at 
an angle of 45° inclined to the Galactic equator from 
(l,b) ps (145°, -10°) to (l,b) (170°, -35°) (Fig- 
ure HI) , consistent with t he boundary of Radio Loop II 
(|Berkhuiisen et all I1973D . A spur of bright polarized 
emission is visible in Figure [5] inside the region with neg- 
ative RM, while strong depolarization is found at the 
boundary where RM changes sign. The Ha image shows 
~ 20° long filaments with a similar orientation with re- 
spect to the Galactic equator, approximately coinciden- 
tal with the lin e where RM changes sign. 

Region B of Sim ard-Normandin fc Kronberg! ([1980), 
centered near I pa 255°, b ss 0°, was identified with the 
Gum nebula by these authors. Structure in RM ampli- 
tude has since bee n noticed in this area in the form of a 
magnetized shell ([Vallee fc Bignell 1 119831: iStil fc Tavlorl 
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Fig. 1. — Figure available as separate jpg file Correlation between RM and Ha intensity from Finkbcincr (2003). Positive RM are 
indicated by red circles, while negative RMs are indicated by blue circles. Circles are scaled according to RM amplitude. 



Fi g. 2. — F i gure available as separate jpg file Correlation between RM and polarized intensity from IWollcbcn ct al. (2006). See 
also Wolleben (2007) for a discussion of the large-scale features in this image. 



2007; that is also visible in Figure [T] We note also that 
there is no clear reversal of the sign of RM across the 
Galactic equat or in the area of the Gum nebula . 

Region C of iSimard-Normandin fc Kronberd ([1980f ) is 
an area with positive RM bounded by 0° < I < 60°, 
0° < b < 60° and was noticed for its wider range in 
RM amplitudes from +50 rad m~ 2 to +300 rad m~ 2 . 
Region C is in the same area as the bottom of Loop I 
(a circle with diameter 116° ± 4° centered around I = 
329° ± 1?5, b = +17? 5 ± 3° Berkhuijsen et al. 1971), but 
ISimard-Normandin fc Kronbergl (|1980l ) considered this a 
coincidence. Figure Q] shows a distinct area inside region 
C with consistently high positive RM in the area 33° < 
/ < 65°, 0° < b < 40°. The median RM in this area 
is 66.7 rad m -2 , while the median RM of the adjacent 
areas 13° < I < 33° and 60° < I < 80°, both with 
0° < b < 40° is 24.6 rad itT 2 . The excess RM amplitude 
over the immediate surroundings is therefore 42 rad m~ 2 . 
The high-latitude boundary of this region of enhanced 
positive RM is not sharp. If it extends well north of 40° 
at the level of ~ 10 rad m~ 2 , it may be responsible for 
the RM excess asso ciated with the Hercules supercluster 
bv lXu et all (pOOfil ). 

Figure Q] also shows enhanced RM amplitude as- 
sociated with some bright H II regions, e.g. at 
(1,6) = (+105°, +4°), (+78°, +2°), (+6°, +23°) (C 
Oph), (-108°, -9°) (Gum), and (-165°, -12°) (A Ori). 
Smaller excesses of RM amplitude are associated with 
more diffuse, low-surface brightness Ha emission, e.g. 
in the Orion-Eridanus region (-180° < I < -135°, 
—55° < b < —5°). We also see a small RM excess asso- 
ciated with the shell of the Orio n-Eridanus superbubble 
(|HeilesHl976tlBrown et all [19951) . 

Figure [3] shows R Ms on an image of 0.2 5 keV diffuse 
X-ray emission from (jSnowden et al.lll997l ). At low lati- 
tudes there is little emission because of strong absorp- 
tion. At high latitudes, emission from the more dis- 
tant halo is observed, while high-latitude clouds absorb 
the halo emission locally, a nd appear as dark silhouette s 
known as X-ray shadows (Snowden et all [19971 119981) . 
Region C' coincides with a large X-ray shadow that has 
an HI counterpart at Vlsr = — 2 km s -1 , and that is part 
of th e outer HI shell of radio loop I ([Berkhuijsen et all 

3. RM STRUCTURE FUNCTIONS 

Structure in RM on angular scales of tens of degrees 
becomes visible because of the high sampling density (> 
1 RM per square degree). The RM density on the sky 
in our data is not sufficient to recognize RM structures 
on angular scales less than a few degrees, but structure 
on small angular scales does affect the variance of RM 
because it introduces differences from one line of sight 
to another. This addition to the variance of RM may 
depend on angular scale, and is conveniently expressed 
as the second-order SF D defined as 



where the sum is over all pairs of sources with a mutual 
distance 59 on the sky, and N is the number of pairs 
included in the sum. The mean distance between an 
RM measurement and its closest neighbor in our data 
is 0?65. SFs can probe RM variance on smaller angu- 
lar scales than this depending on the number of source 
pairs available to evaluate the variance with reasonable 
confidence. For our analysis, the smallest angular scale 
is taken to be the scale for which we have a minimum 
of 20 source pairs, or 0?1, whichever comes first. The 
largest angular scale we consider for the SFs is 59 = 10°, 
so the RM gradient related to the large-scale magnetic 
field does not contribute significantly to D. 

All possible sources of variance contribute to the am- 
plitude of the SF. We may write D as a sum of variance 
terms that we consider statistically independent 

D = <?L + <Tl 2 GM («) + ^ISmW + ^oise(^), (3) 

with af nt the variance of RM originating inside the 
source, of GM the variance in RM originating in the in- 
tergalactic medium, for different lines of sight separated 
by angle 59, er 2 SM the variance of RM originating in the 
interstellar medium for lines of sight separated by angle 
59, and cr 2 oisc the variance resulting from measurement 
errors and noise. 

The variance c 2 oise should be subtracted from D to ob- 
tain the sum of the three te rms of astrophysical interest 
(e.g. lHaverkorn et all 120 04). This term was calculated 
from the errors on the RMs contributing to pairs sep- 
arated by angular distance 59. Each term in the sum- 
mation of Equation [5] contains differences RM m — RM n 
(m > n indices of data in the RM catalog). For each 
pair of sources included in the summation, the variance 
added by noise is cr 2 ^ + c 2 where a m and a n are the 
standard errors for RM m and RM n respectively. This 
assumes that the error bars in the RM data base provide 
a realistic estimate of the measurement errors. We defer 
this discussion to Section [4] In principle c 2 oisG can de- 
pend on 59, but for the SFs presented here, it is virtually 
constant. The median error in our complete RM catalog 
is ~ 10.8 rad m -2 , so the error term can be a signifi- 
cant source of variance and it affects the slope of the SF 
on small angular scales. Unless stated otherwise, SFs in 
this paper or derived properties such as amplitude and 
slope have been corrected for variance associated with 
measurement errors. 

Figure O shows a simulation of the effect of noise on 
a SF, and the noise-subtracted SF. The simulated data 
consisted of 1978 points drawn at random from an image 
with a random field with input SF with slope 0.6 kindly 
provided by V. Uritsky. For Figure [S^,, measurement er- 
rors for each data point were simulated by first drawing a 
standard deviation from a uniform distribution between 
0.4 and 2.0 (data units), then drawing a number from 
a Gaussian distribution with the chosen standard devi- 
ation, to be added to the actual data value. The stan- 
dard deviation itself was recorded as the error bar for 
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Fig. 3. — Figure available as separate jpg file Correlation between RM and 0.25 keV diffuse X-ray intensity from Snowdcn ct al. 
(fT998h . 



Fig. 4. — Figure available as separate jpg file RM structure on the sky centered at (I, b) = (180°, 0°), leaving out data with 
\RM | < 25 rad m — 2 , with positive RM in red and negative RM in blue, and no scaling of the circles according to RM amplitude. This 
representation is analogous to showing a the lowest contour in a radio map at the level of a few sigma. It visualizes RM structure and 
boundaries where the sign of RM changes. Regions with a high RM amplitude appear more densely sampled. Note the arc-like structures 
in RM on angular scales of tens of degrees that add a few hundred rad 2 m -4 to the SF amplitude on large angular scales. 
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Fig. 5. — Simulation of the effect of noise on a SF, assuming two 
different levels of RM errors. The upper solid curves with dots rep- 
resent the measured SFs with noise. The dashed curves show the 
power calculated from the adopted errors. The points with error 
bars represent the SF after subtracting the noise power. The SF 
made with the original 1978 noiseless data points is represented by 
the lower solid curve. The measurement errors in (b) are on aver- 
age two times as large as the measurement errors in (a), resulting 
in a noise power that is 4 times higher in (b). The relative strength 
of the noise power in (a) is more representative for the data in this 
paper. 

each RM. The figure shows the SF as observed (upper 
solid curve with dots), the noise power cr 2 oise evaluated 
from the error bars of RM (dashed curve), and the SF 
corrected for noise power (points with error bars), com- 
pared with the SF of the same 1978 noiseless data points 
(lower solid curve). For Figure [SJd, the same set of origi- 
nal data points was used, but the standard deviations of 
the RM errors were drawn independently from a uniform 
distribution between 0.8 to 4.0 data units. This on aver- 
age doubles the measurement errors compared with the 
situation shown in Figure [S^,, and quadruples the level of 
the noise power in Equation [3] 

We retrieve the input SF in the simulation shown in 
Figure [S] and other simulations that included SFs of vari- 
ous slopes and noise levels. These simulations confirmed 
our evaluation of the noise power even though the errors 
in individual RMs varied by a factor 5. On the smallest 
angular scales, we see occasional large variation in the 
data after noise subtraction. This is to be expected, be- 
cause the contribution of the noise to the total variance 
is large, but our estimate of the noise power is uncertain 
because relatively few source pairs with small angular 
separation exist in the catalog. We can easily identify 
such points and remove them from consideration. 

The subtraction of the SF of the noise is a statistical 
correction that is more accurate for larger samples. The 
SF of the measurement errors is not necessarily flat, for 
instance in the case of periodic variations in the noise of 
a mosaicking survey with widely separated field centers, 
or in the case of a deep survey of a small area of the 
sky where the noise increases towards the edge of the 
survey area (affecting all source pairs with large angular 
separation). In our analysis we do not assume that the 
SF of the noise is flat, or that all RM error are equal. We 



calculate the SF of the noise for each 58 from the actual 
error bars of all included RMs. However, we found no 
systematic variati on of the SF o f the noise with angular 
separation for the iTavlor et al.l {2009) RM catalog. 

3.1. Direction dependence of the RM SF 

Before we derive an all-sky image of SF parameters, 
we examine SFs for a few areas of interest inspired by 
results discussed in the previous section. These regions 
visualize the significance of differences in SF parameters 
with direction. Table [T] lists these regions, and SFs are 
shown in Figures [6j [71 and [8] 

Experience from t h e literatu r e (|Haverkornl 120071 : 
lHaverkorn et al.1 120081: IRov et al.l l2008t l and from the 
present data (|Still I2009T) shows that the slope of the SF 
is not always constant, but sometimes changes slope at 
an angular scale 58 w 1°. The precise angular scale at 
which a break in the power law occurs is not necessarily 
the same everywhere, but visual inspection of SFs across 
the sky suggests that if a break occurs, the SF slope is 
considerably higher for 58 < 1° than for 58 > 1°. We 
therefore fit separate power laws to the SFs at angular 
scales 58 < 1° and 58 > 1°: 

D{58) = A, 68 a * (4) 

with i = 1 for 56 < 1° and i = 2 for 58 > 1°. Table CD 
lists results of such fits for the SFs shown in Figures [51 [71 
and [U The fits for small angular scales are much less 
constrained because of the limited number of close pairs, 
and a relatively larger contribution of cr^ oiso to the total 
variance on small angular scales. We use the fit param- 
eter A2 as a measure of the amplitude of the SF. A2 is 
the variance of RM for an angular scale of 1°, corrected 
for variance due to measurement errors. 

3.1.1. Polar regions 

Figure [5] shows RM SFs within 15 degrees of the North 
Galactic Pole (NGP; 1019 sources) and the South Galac- 
tic Pole (SGP; 752 sources). Parameters for these SFs are 
listed in Table [U We have fewer RMs in the SGP region 
because of the declination limit of the NVSS (5 > -40°), 
and because of less sky coverage with boths IFs required 
to derive a RM. The contribution of cr^ oisc is indicated 
by the dashed curves, and subtracted in Figure [Hb and 
d. Subtracting the noise variance that is approximately 
independent of angular scale 58, reduces the amplitude 
of the SF significantly, and also steepens the slope of the 
SF on small angular scales. 

The amplitudes of the SFs at the NGP and the SGP 
are signifcantly different. The amplitude of the SF in Fig- 
ure [6^, on angular scales more than 2° i s consistent with 
the am plitude ~ 480 rad 2 m~ 4 found bv lSimonetti. et all 
lfl98l for this region. These authors did not subtract the 
CT noise t erm in Equation[3] The SF at the NGP also shows 
a slight increase in amplitude from angular scales 0?2 to 
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TABLE 1 

Regions and parameters for SFs shown in Figures [6] [7] and [8] 



Region 


imin 'max 6 m in Vax log(j42) 


on 


a 2 


Figure 




(degr.) (degr.) (degr.) (degr.) 









NGP 


-180 


180 


75 


90 


SGP 


-180 


180 


-90 


-75 


Orion 


-180 


-135 


-50 


-10 


anti-Orion 


-180 


-135 


10 


50 


Region C' 


33 


68 


10 


35 


anti-Region C' 


33 


68 


-35 


-10 


Region A' 


80 


150 


-40 


-20 


anti-Region A' 


80 


150 


20 


10 


Fan 


115 


160 


5 


25 


anti-Fan 


115 


160 


-25 


-5 


Gum 


100 


100 


20 


20 



2.40 (0.03) 
2.67 (0.03) 
2.92 (0.02) 

2.67 (0.01) 
3.06 (0.01) 
3.15 (0.03) 
2.81 (0.01) 

2.68 (0.01) 
2.78 (0.02) 
2.98 (0.02) 
3.72 (0.04) 



0.63 (0.39) 
0.99 (0.58) 
0.79 (0.08) 
1.10 (0.39) 
0.33 (0.20) 
0.31 (0.48) 
0.82 (0.30) 
0.74 (0.26) 
-0.02 (0.12) 
0.77 (0.17) 
1.08 (0.41) 



0.05 (0.04) 
0.02 (0.04) 
0.50 (0.03) 
0.31 (0.01) 
0.39 (0.02) 
0.58 (0.05) 
0.51 (0.02) 
0.38 (0.03) 
0.58 (0.03) 
0.59 (0.04) 
1.14 (0.09) 
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Fig. 6. — (a) RM SF in the NGP region (b > 75°). The dashed line gives the power in the SF from measurement errors, (b) SF for the 
NGP region with noise power subtracted, (c) Same as (a) for the SGP region (b < —75°). (d) Same as (b) for the SGP region. 



10°, but the error of a-i in this region is large (Table [T]). 
iMao et al. I (|2010l ) reported a slope 0.08 ± 0.01 in the 
NGP region and 0.03 ± 0.01 in the SGP region. Our SF 
at the SGP appears flatter than at the NGP. 

A difference in the slope and amplitude of the SFs 
of the two Galactic pole regions suggests a significant 
Galactic foreground in these regions. Assuming of GM 
is negligible, and of nt is independent of 89 for unre- 
lated sources, a rise in the noise-corrected SF or dif- 



ferences in amplitude must be related to differences in 
a iSM- Since the amplitude of the SF at the SGP is 
higher than at the NGP, and the SF at the NGP ap- 
pears to have a positive slope on smaller angular scales, 
it follows that both Galactic poles have significant fore- 
ground RM variations. An upper limit eri nt < 10 rad m -2 
is suggested by the RM variance in the NGP region on 
angular scales less than 1°. We do not expect a- lnt to 
be much less than 5 rad m -2 because the variation of 
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Fig. 7. — (a) RM SF in the Orion region (black) and its northern counter part ("anti-Orion"; grey) as defined in Table [T] The noise 
variance crj? . has been subtracted, (b) SF for Region C' (black) and its southern counterpart "anti-region C" (grey), (c) SF for Region 
A' (black) and its northern counterpart "anti-region A'" (grey), (d) SF for the Fan region (black) and its southern counterpart "anti-Fan" 
(grey). 



RM across individ ual sources a t high latitude on angular 
scales less than 1' (|Leahvlll987D is not much smaller than 
our upper limit determined at angular scales 10' - 25' 
from Figure Eb, while mod els of the foreground suggest 
lower values ([Sun fc Reich! [2001 ) . 
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additional variance arises in a single Faraday screen, the 
rms RM amplitude of this screen would be 14.7 rad m~ 2 , 
but there is no evidence that the higher variance results 
from a separ ate additional Farad ay sc reen in the southern 
hemisphere. iTavlor et ail (|2009[ ) and lMao et al. I (|20100 
reported that the mean RM near near the SGP is larger 
than near the NGP. The difference in RM variance found 
here is consistent with the difference in mean RM in the 



sense that both variance and mean are larger by a factor 
~ 2 in the south. 



3.1.2. RM SFs at lower latitude 

Figure [7] and Figure H] show SFs of selected 
intermediate-latitude areas defined in Table [TJ The pairs 
of SFs in each panel of Figure [7] represent regions with 
the same longitude, mirrored with respect to the Galactic 
plane. The Gum region is shown separately in Figure [S] 
because its mirror region is not covered by the NVSS, 
and the amplitude of the SF is much larger than the SFs 
shown in Figure[7J The first impression of Figure[7Jis that 
the amplitude of the SFs is not symmetric with respect 
to the Galactic plane. The highest SF in each panel of 
Figure [7] is in the southern Galactic hemisphere, showing 
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the same asymmetry as Figure [5] On average these SFs 
suggest that the variance of RM in different areas in the 
southern Galactic hemisphere is ~ 50% larger than in 
the northern Galactic hemisphere. iTavlor et al.l (j2009f ) 
found that the mean RM across the southern Galactic 
hemisphere is consistently higher than in the northern 
Galactic hemisphere. We will show that the asymmetry 
in RM variance is also global. 

Figure [7] and Figure [3] also show significant differences 
in the slope of the SFs. The SF of the Gum region is 
among the steepest identified in the data. The Orion re- 
gion has a significantly steeper slope than its counterpart 
region across the Galactic plane, and the same is true for 
Region A'. In other areas, such as the Fan region, we see 
a difference in amplitude but no significant difference in 
slope. 

Figure|9]the sky distribution of A2 and oli derived from 
a dense grid of SFs, using sources inside a circular aper- 
ture with diameter 15°. The slope and amplitude of the 
SFs for 69 > 1° have a similar accuracy as those listed 
in Table [TJ because of the similar number of sources in- 
volved. The uncertainty in the amplitude is therefore 
approximately 5%, and the uncertainty in the slope is 
approximately 0.03. Most of the structure visible in Fig- 
ure [5] represents significant differences in amplitude or 
slope of the RM SF. 

The asymmetry in SF amplitude between the NGP and 
SGP shown in Figure |6] is also visible in the top panel of 
Figure [H There are only a few areas in the south where 
the amplitude drops below 300 rad 2 m~ 4 , while most of 
the NGP region has amplitudes less than 300 rad 2 m~ 4 . 
Towards the Galactic plane, RM fluctuations increase. 
The largest fluctuations (4.8 x 10 5 rad 2 m -4 ) are found 
in the Galactic plane, in the direction I f» 90°. RM fluc- 
tuations may be approximately symmetric with respect 
to I = 0° , but the declination limit of the NVSS removes 
the fourth Galactic quadrant and part of the third Galac- 
tic quadrant. 

The amplitude of the SFs in Figure [9^ traces the 



brightest extended Ha emission, with significantly higher 
amplitudes towards the inner Galaxy than towards the 
outer Galaxy. This contrast shows the excentric loca- 
tion of the Sun in the Galaxy, and thus suggests that 
the low-latitude RM variance on an angular scale of 1° 
originates from a line of sight distance of several kpc. 
Structure at higher latitudes are probably more local in 
nature. Around longitude I — —110°, the contours of SF 
amplitude follow the outline of the Gum nebula in the 
Ha image. The low-latitude peaks follow bright emission 
in the Galactic plane, except perhaps the local maximum 
near I = 40°. The brighter parts of the Orion- Eridanus 
region (I = —170°, b = —40°) are also associated with 
enhan ced RM variance. Th e outlines of radio Loops I- 
IV (jBerkhuiisen et al.|[l973h are also shown in Figure |H 
Only Loop III may have have a counterpart in SF am- 
plitude. 

The slope of the SF (Figure [9)d) is highest in the Galac- 
tic plane, but it does not show the same concentration 
toward the Galactic plane as the amplitude. The largest 
slopes are a?, rj 1.3 in the direction of the £Oph pjn 
region (I = 5°, b — +24°), and the Gum nebula. At 
higher northern Galactic latitude, we see a systematic 
trend that the slope of the SF is shallower towards the 
Galactic center than towards the Galactic anticenter at 
the same latitude. The contour of = 0.2 reaches 
b ~ 50° towards 1 = 0°, while it is found around b « 80° 
towards I = 180°. It is not clear whether the same behav- 
ior also occurs at southern Galactic latitudes, because of 
the unsampled region S < —40°. The a2 = 0.2 contour is 
found mostly between b — —60° and b = —70°, except for 
0° < I < 30° where it is closer to the Galactic plane simi- 
lar to its northern counterpart. The over- all distribution 
of the SF slope at high latitudes is more symmetric with 
respect to the Galactic plane than the amplitude. There 
is a suggestion of a steeper slope in the rim of radio Loop 
III. 

Figure [TU] shows the amplitude and slope of RM SFs 
as contours on po larized intensity of di ffuse emission 
at 1.4 GHz from iWolleben et"afl (|2006l ). The main 
features in this polarized intensity were described by 
IWolleben I (|2007[ ). The gradual trend in slope if high- 
latitude SFs with longitude mimics the bound aries of the 
high-l atitude polarized emission described by IWolleben I 

pool . 

In the first Galactic quadrant (0° < I < 90°), the drop 
in polarized intensity at b w 30° and b ~ —30° noted 
bv IWolleben et afl (|2006l ) corresponds with RM variance 
of 600 ±100 rad 2 m~ 4 north and south of the Galactic 
plane, although we see extended polarized emission in 
some areas with RM variance up to 900 rad 2 m~ 4 . The 
polarized intensity drops abruptly by a factor ~ 3 at the 
boundary, while the average percentage polarization is a 
factor ~ 10 smal ler at lower latitudes than for b > 30° 
IjWolleben 1 120071 ) . Some of this difference may result 
from depolarization at a la rger distance. T he 23 GHz 
WMAP polarization image (jPage et al.l 120071 ) shows the 
North Polar Spur polarized at lower latitudes, implying 
some depolarization occurs for \b\ < 30° at 1.4 GHz. We 
take the abrupt decrease of a factor ~ 3 in polarized 
intensity as a measure for the amount of depolarization 
by a possible foreground screen. 

In the second Galactic quadrant (90° < I < 180°), 
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Fig. 9. — Figure available as separate gif file Direction dependence of SF amplitude A2 and slope a-z- Top: amplitude for 86 = 1° . 
Contour levels are 300, 600, 1200, 2400, 4800, 9600, and 19200 ra d 2 m~ 4 . Th e lowes t contours show tickmarks in the direction of smaller 
amplitude. Grayscales represent He? intensity in Rayleigh from Finkbciner (2003) on a logarithmic scale indicated by the color bar. 
The closed black l oop m arks the lower southern boundary of the NVSS (5 = 40°), and the gray loops indicate radio loops I- IV from 
IBerkhu iiscn ct al. (1973). Bottom: Slope of the SF on scales > 1° («2). Contour levels range from 0.2 to 1.2 with increments 0.2. Contours 
at the levels 0.2, 0.4, and 0.8 are labeled. The grayscales are the same as in the upper panel. The SF slope varies between —0.07 and 1.28. 



Fig. 10. — Fi gure available as sep arate gif file Same as Figure[9] with gray scales representing polarized intensity of diffuse emission 
at 1.4 GHz from Wollcbc netaLl 12009) . 



the brightest polarized emission is observed at lower lat- 
itudes north of the Galactic equator in the Fan region. 
The brightest polarized emission of the Fan region is in- 
tersected by the 900 rad 2 m -4 RM variance contour. In 
contrast to the first quadrant we do not see a clear anti- 
correlation between polarized intensity and RM variance, 
with the possible exception of Radio Loop III. 

The third Galactic quadrant (180° < I < 270°) is only 
partially sampled by the present data. The region of 
the Gum nebula as traced by the 600 rad 2 m~ 4 con- 
tour appears depolarized. A narrow polarized feature 
between longitudes —180° > I > —155°, and latitude 
— 10° < b < 0° is an extension of polarized emission in 
the second quadrant. It does not appear to be correlated 
with structure in RM variance. 

4. DISCUSSION 

4.1. RM contributions from the local ISM 

RMs from polarized extragalactic sources probe the en- 
tire line of sight from the source to the observer. Most 
extragalactic sources experience relatively little internal 
Faraday rotation, and intergalactic Faraday rotation is 
believed to be substantially smaller than the intrinsic 
RM variance of the extragalactic sources. With the ex- 
ception of some sources with high internal Faraday rota- 
tion, most of the RM originates from the section of the 
line of sight inside our own Galaxy. RM amplitude also 
increases significantly towards the Galactic plane. Al- 
though the RM is an integral along the line of sight, a 
single object can contrib ute significantly. Th e Gum neb- 
ula (distance ~ 400 pc iBrandt et al.1 119711 ) contributes 
significantly to the RM, even though the line of sight 
towards the Gum nebula intersects the inner Galaxy, 
where the elect ron density is high o ver a line of sight 
of several kpc (jCordes fc Laziol [20021) . The Orion star 
for ming region, also at a distance of approximately 400 
pc (|De Zeeuw et al.l Tl999f ). is at higher latitude, so its 
background is a relatively short line of sight through the 
Galaxy. 

The line of sight through the Galaxy is much shorter 
for high-latitude sources, and structures in the local ISM 
may affect the RM significantly. Some high-latitude H 11 
regions were identified as such, consistent with earlier 
results that showed that ind ividual H II region s can have 
a significant impact on RM ( Mit ra et al.l l2003) . We have 
found evidence that Region C is associated with a nearby 
cloud revealed by its absorption of soft X-rays from the 
Galactic halo that is at least partially neutral because it 
is visible in H 1 emission. This associates region C' with 
the outer HI shell of Loop I. The enhanced RM amplitude 
in region C' spans only 35° in Galactic long itude, with its 
centr al longitude around I = 45° (see also lStil fe Tavlorl 
l2007lL 

Towards the Galactic anticenter, where we expect al- 



most no contribution from the large-scale Galactic mag- 
netic field, we find curved filamentary structures in RM 
with amplitudes of a few tens rad m~ 2 , positive and nega- 
tive (Figure QJ. From Equation [1] we obtain an order-of- 
magnitude estimate for the RM imposed by a structure 
along the line of sight with electron density n e , internal 
regular magnetic field B, and line-of-sight extent I: 

RM « 8 f Ue , ) ( ( —^-) rad m" 2 . (5) 
V0.1 cm- 3 / VI ^gJ VIOOpc/ w 

The RM amplitudes and morphology of these structures 
resemble one or a superposition of several bubble walls 
with size of the order of ^100 pc (with a large uncer- 
tainty). One of these is the Orion-E ridanus superbubbl e 
(170° < I < 220°, -30° < b < -50° (|Brown et al.lfl995l) . 
visible in emission in Figure [T] and in absorption in Fig- 
ure 02 These filaments contribute to the amplitude of 
RM SFs on large angular scales, increasing the SF slope 

The overarching picture that emerges is one of RM 
structure in the local (< few kpc) ISM on large scales is 
mostly associated with bubble walls, H 11 regions, and lo- 
cal clouds, region A' (Tabled]) being a possible exception. 
If this is representative for an average few kpc line of sight 
through the Galaxy, the RM variance added by struc- 
tures on scales of a few hundred pc should contribute 
significantly to the total variance in RM at low Galactic 
latitudes. If bubble walls contribute significantly, mod- 
els of RM variance should include a correlation between 
field strength and density on scales of a few hundred pc. 

4.2. Amplitude and slope of SFs 
4.2.1. Comparison with Mao et al. (2010) 

IMao et al. I (pOlOh derived SFs at the NGP and SGP 
caps (latitu de \b\ > 77 ° ) bas ed on 354 and 319 RMs 
respectively. IMao et al. I (|2010f ) find almost identical SF 
amplitudes at both Galactic poles, an apparent contra- 
diction of the result in Figure [6] These authors also 
presented a comparison of their RMs with our RM cat- 
alog, and foud a small correlation coefficient between 
the two samples (0.39 towar d s the NGP, and 0.36 to- 
wards the SGP). IMao et al. I (|2010l ) suggested that the 
differences between the samples are related to multiple 
RM components in some sourc es that canno t be re solved 
in the two- frequency data of I Taylor et alj ((2009), but 
they did not provide specific examples where this would 
be the case. While it is true that the two-frequency 
data in I Taylor et al.l (|2009T) cannot resolve non-linearity 
in the relation between polarization angle and A 2 , ro- 
tation measure s ynthesis with limi ted A 2 coverage, as 
in the sample of IMao et al. I (I2010D has its l imitations 
(|Brentiens fc De Bruvnll2005t iFrick et al. II2010D . It was 
recently pointed out by L. Rudnick et al. (unpublished) 
that RM synthesis can provide erroneous solutions in a 
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Fig. 11.— Comparison of RMs from ITavlor et al.l | |200S|) with RMs from IMao et al. I (l2010h . divided by stated RM error, (a): NGP 
region with RM error < 8 rad m — 2 in both samples, (b): NGP region with RM error > 8 rad m -2 in ITavlor et all ||2009| 1 (x), and RM 
error > 8 rad m — 2 in Mao ct al!] 120101) (+). If the errors in both surveys are > 8 rad m — 2 both symbols appear, resulting in a *. (c): 
same as (a) for the SGP region, (d): same as (b) for the SGP region. 

ter in Figure [TTJ but still raise the amplitude of SFs de- 
rived from both data sets. 

Figure [11] a and c both show significant RM signal. 
The scatter in these relations provides an estimate of the 
variance introduced by measurement errors. We examine 
the subsets with small errors shown in Figure [TT] a and 
c, as well as the complete sample. For the subsets with 



simple system with just two discrete RM components. 

In order to gain a better un derstanding of the dif- 
ference between SFs derived by IMao et al. I (|2010ft and 
those derived here, we must first compare the RM sam- 
ples in more detail. The median error for sources com- 
mon t o both RM catalogs is 8 rad m -2 in Z 
(2009), 5 rad m" 2 in the NGP region of 
2 in the SGP region of 




and 3 rad m 



'avlor ct al 



Mao et al. 



Mao et al. 



small errors, the total variance is 90 rad m in the 



Figure [UJ shows scatter plots of RMs matched north and 135 rad 2 m 4 in the south. The error bars 



between the two samples, subdivided into RMs with er- 
ror < 8 rad m~ 2 (panels a and c), and all other RMs 
(panels b and d). 

When considering only RMs with small errors, we see 
a correlation between the two datasets at both Galac- 
tic poles (Figure fT T] a and c) . The la rger scatter in the 
SGP area found bv lMao et all (120101) is relat e d to R Ms 
with larger error bars, in either [Taylor et all (|2009f ) . or 
IMao et al. I ([20101 ) , or both. We do not calculate a corre- 
lation coefficient, because the correlation coefficient be- 
tween independent RM surveys of low-RM sources must 
be small. Unresolved multiple RM components may af- 
fect both surveys in a similar way, and thus introduce 
correlated systematic errors without increasing the scat- 



m the 
4 and 



(added in quadrature for each source) account for a total 
variance of 37 rad 2 m -4 in FigurefTTa and 36 rad 2 m~ 4 in 
Figure [Tib. For the complete sample, the total variance 
is 151 rad 2 m -4 in the north and 193 rad 2 m 
south, while the error bars account for 103 rad 2 
90 rad 2 m~ 4 respectively. 

The stated errors underestimate the true variance in 
the correlation of RMs common to both samples, but 
more so if the state d errors are small. The RM errors in 
ITavlor et all ([200l were estimated by error propagation 
of the n oise in the Stokes Q and U images, while the 
errors in IMao et al. I ((2010) were estimated as half the 
width of the RM spread function divided by the signal 
to noise ratio in polarization. The error estimates in both 
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surveys converge to zero as the signal to noise ratio in 
polarization increases without bound. Perhaps the true 
errors do not converge as rapidly as assumed. Before we 
consider the effect of this underestimation of errors on the 
amplitude of SF at the Galactic poles, we first consider 
another aspect of the error statistics of the data. 

Figure [T2] shows the stated errors as a function of 
RM amplit u de fo r the NGP and SGP regions of the 
iMao et al. I (|2010[) sample. Black dots represent RMs 
used in the construction of SFs, and gray crosses resprc- 
sent rejected RMs. While the median error in Figure [T"2"a 
is larger, Figure \V2b contains data with large RM errors 
that contribute to the increased scatter in Figure [TTH . 
RMs with amplitude more than 20 rad m~ 2 in Figurc[T2b 
are much more likely to have a large error, introducing 
a correlatio n between RM erro r and RM amplitude: the 
SGP data of lMao et aT~l ([2010D are heteroskedastic. The 
highest RM amplitudes (a few percent of the total sam- 
ple) are rejected in both regions, but the distribution of 
rejected RMs is more skewed toward high RM amplitude 
in the south (F igure [12b ). 

Some RMs in Ma o et al. I ()2010l ) were rejected in order 
to eliminate AGNs with strong internal Faraday rotation, 
and to minimize the effect of known Galactic and ex- 
tragalactic structures. This approach is preferred if one 
wishes to determine the intrinsic RM variance of AGN- 
powered radio sources, or investigate small-scale turbu- 
lence at high Galactic latitude. This paper deals with the 
complete RM foregound. Figure IT51 shows the RM errors 
in relation to RM amplitude for our data in the NGP and 
SGP regions. The larger RM errors in our data are com- 
pensated by the ~ 3 times larger sample size, especially 
when probing fluctuations on small angular scales. 

What does this all mean for the SF amplitudes at the 
Galactic poles? Figure[Hshows RM SFs at the NGP and 
SGP. The g r ay SF s were made with the censored data of 
IMao et al. I (|2010D (354 RMs at the NGP and 319 RMs at 
the SG P). We find a smaller amplitude than lMao et al. I 
(2010), because these authors underestimated the noise 
power by approximately a factor 2. The problem was 
identified and the authors confirmed they recover the 
SF shown in Figure [14] (A. Mao, private communica- 
tion). The black points in F igure IT41 sh o w SF s made 
with the complete sample of IMao et al. I ((20 10 s ) . The 
origin of the large scatter for 59 < 1° in Figure [14a was 
discussed in Section Ep (Figu re El). Using the complete 
sample of M ao et al. I ( 20101 ) , we confirm that the SF 
in the SGP region is a factor ~ 2 higher than in the 
NGP region, but both SFs in Figure Q3] have a some- 
what smaller amplitude than those in Figure [SJ Both 
data sets suggest a small upward slope of the SF in the 
NGP region In T able CD we list a 2 = 0.05 ± 0.04, while 
IMao et al. I (|2010t ) found 0.08 ± 0.01. For the SGP region 
we find a 2 = 0.02 ± 0.04, while IMao et~al~l (poToT ) found 
0.03 ± 0.01. 

If we assume that the noise power in our SFs is really 
a factor 1.5 higher than that calculated from the RM 
errors, we would retrieve the black curves in Figure 1141 
This would be equivalent with the RM errors bein g un- 
derestimated by a factor 1.22 in lTaylor et al.l ()2009l ). and 
account for the variance in the comparison of RMs from 
the two surveys. While this appears a plausible expla- 
nation of the difference, the actual comparison may be 
more complicated. This analysis assumes that Figure [TT] 



reveals all RM errors, while it does not reveal possible 
systematic errors common to both surveys that may in- 
crease SF amp litude. Another qu estion is the nature of 
the sources in IMao et al. I (|2010f ) for which no reliable 
RM could be found (6% of their data in the NGP region 
and 28% in the SGP region. 

It is encouraging that we find consistency between SFs 
from independent data sets, after accounting for all noise 
power traced by direct comparison of RMs common to 
both surveys. This consistency requires just a 22% in- 
crease in the errors quoted by Taylor et al. (2009|) (in 
a low-RM environment), corresponding with only a few 
rad m -2 . It suggests that eit her the effect of mu ltiple RM 
components on the da ta of ITavlor et aL (|2009t ) is quite 
small, or the data from lMao et al. I (|201C ) are affected in 
the same way for most sources. However, such a system- 
atic error should still increase the amplitude of the SFs, 
assuming it is uncorrelated between sources. From the 
total amplitude of the SFs at high latitude, we estimate 
that any such systematic error should be smaller than 
~ 10 rad m~ 2 . 

Our analysis shows that both RM surveys point indi- 
cate a larger RM variance in the direction of the SGP 
than in the direction of the NGP. Our SFs at inter- 
mediate Galactic latitude show the same trend, and 
these are less sensitive to the amplitude of the noise 
power beca use it contributes mu ch le ss to the tot a l vari - 
ance. Both ITavlor et all ([20091 ) and IMao et al. I (|2010l ) 
reporte d a higher mean R M in the south than in the 
north. ITavlor et al.l (|2009l ) found this to be the case in 
the entire southern Galactic hemisphere. 

4.2.2. Distribution of RM variance over the sky 

Variance in RM originates from fluctuations in electron 
density and magnetic field on a range of scales in the ISM, 
including H n regions, bubbles, turbulence, and possibly 
also molecular clouds with a small ionized fraction but 
a strong magnetic field. The variance in RM should be 
approximately proportion al to the length of t he line of 
sight through the Galaxy (jSokoloff et al.lll998[ ). 

Figure[9]shows that the amplitude of SFs on an angular 
scale 68 — 1° follows the distribution of Ha intensity in 
the Galactic plane. The highest amplitudes at low Galac- 
tic latitude are found in the direction of the Local Arm 
(/ ?» 90°) and other regions along the Galactic plane. 
The peak in RM variance at (l,b) ~ (40°, 0°) occurs in 
a region of the Galactic plane with strong extinction at 
visible wavelengths. The enhance variance in the Gum 
region is clearly extented in Figure [S] with morphology 
similar to Ha intensity. The low-latitude RM variance is 
smaller in the direction of the Galactic anti-center, where 
the line of sight through the Galaxy is shorter. These ob- 
servations are consistent with the expectation that vari- 
ance in RM increases with the length of t he line of sight 
through the Galaxy (jSokoloff et al.l ll998). Higher elec- 
tron density and stronger turbulence related to more in- 
tense star formation in the inner Galaxy probably also 
contribute to the high RM variance for lines of sight 
through the inner Galaxy. The outer scale in RM vari- 
ance associated w ith inter arm regio ns, b ut not with spiral 
arms r eported bv lHaverkornl (|2007l ) and lHaverkorn et "aTl 
(2008) indicates that besides line-of-sight distance, the 
local conditions are also important. 

The difference in RM variance between the Galactic 



Structure in the Rotation Measure Sky 



11 



30 



20 



00 



1 o ■ 



30 



20 



(b) 



20 40 60 80 100 
|RM| (rod m" 2 ) 




30 



20 40 60 80 100 
|RM| (rod m" 2 ) 



Fig. 12, — Correlat i on of RM error with RM amplitude for the 
sample of Mao et al. (201J5). (a): NGP region observed with the 
WSRT, (b): SGP obser ved with the ATCA . Black dots represent 
RMs used in the SF by IMao et al71 (120101 ) (heavily saturated by 
overlapping symbols in so me parts of the figu re). The gray crosses 
represent data rejected by IMao et al. 1 {20jJJ) for use in SFs. 

northern and southern hemispheres (Figure |H] and Ta- 
ble [T]) is particularly interesting. Assuming the intrin- 
sic RM variance of extragalactic sources to be a- mt = 
7 rad m -2 (Equation [3]), the RM variance introduced by 
the Galactic foreground is 202 rad 2 m~ 4 at the NGP, 
and 419 rad 2 m" 4 at the SGP (120 rad 2 m" 4 and 307 
rad 2 m~ 4 respectively using SFs from Figure [T4"|) . The 
RM variance near the SGP originating from the Galactic 
foreground is therefore a factor 2 larger than the Galac- 
tic RM variance near the NGP. The difference in vari- 
ance is equivalent with an additional source of RM varia- 
tion with amplitude 14 rad m -2 in the southern Galactic 
hemisphere. The Sun is located 20 pc from the Galactic 
mid-plane in the direction of the NGP. This distance is 
too small compared with the scale height of the ionized 
gas to make a significant difference in RM variance. It 
is more likely that we see a local asymmetry of density 
or ma gnetic field, o r both with respect to the Galactic 
plane. TTavlor et al.l (|2009D also found a factor ~ 2 dif- 
ference in the mea n RM of the NGP and the SGP, while 
IMao et aT~l (|2010l ) found a significant positive RM at the 
SGP but not at the NGP. 

The slope of the SF is higher near the Galactic plane 
with most values in the range from 0.4 to 0.8 and ex- 
cesses up to 1.3 in the direction of two nearby Hn re- 
gions (C Oph and the Gum nebula). The slope of our 
SFs (Table [T]) is generally smaller th an the 5/3 slope 
expected for Kolmogoroy turb ulence (|Armstrong et all 
Il995t IMinter fc Spanglerlll996f) . to t he sm a llest a ngular 
scales probed by our data (~ 0?3). IBeckl (|2007l) found 
RM SF slopes ~ 0.3 in the nearly face-on galaxy NGC 
6946 on scales up to 6 kpc, which is in the range of values 
of a 2 in Table [TJ 

If a change in slope occurs in the range of angular scales 
probed by the current data (0?1 < 59 < 10°), the slope 
on larger angular scales is usually smaller. A change in 
SF slope is sometimes interpreted as a physical scale in 
the turbulence probed by RM variance, that could be 
an outer scale length if the SF becomes fiat on large 
scales. An implicit assumption is that the RM variance 
traces stochastic variations of electron density and mag- 
netic field related to turbulence in the magneto-ionized 
ISM. If deterministic structure, e.g. an H II region that 
may or may not be detectable in current Ha surveys, or 
a meso-scale structure in the magnetic field, contributes 
significantly to the variation of RM with position over the 
area for which the SF is evaluated, the shape of the SF 
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Fig . 13. — Same as Figure 1121 for the data from I'l'avlor et al.l 
(2009). (a) NGP region, (b) SGP region. 

may reflect the angular scale of the ionized structure, not 
the turbulence along the line of sight. The present data 
suggest that this may be true in some directions. The 
steepest SFs at low latitude are associated with bright 
extended H n regions in Figure [51 Figure @] also shows 
RM structure on angular scales > 10° with amplitude of 
a few tens of rad m -2 , sufficient to affect the slope of 
the SF between 59 = 1° and 69 = 10°. Sampling with 
a denser RM grid can reveal such structures, as in th e 
case of the RM anomaly in Cygnus (jWhiting et al.l[2009l) . 
Future wide-area RM surveys with the Arecibo radio 
telescope (GALGACTS), the Australian SKA Pathfinder 
ASKAP (POSSUM), and with the Square Kilometre Ar- 
ray will provide denser sampling that may reveal under- 
lying structure that affects the slope of our SFs, while 
not recognizable i n the presen t data . 

Model SFs by ISun et all (12001) and ISun fc Reichl 
(|2009f) that assume a turbulent Kolmogorov power spec- 
trum for density and magnetic field, predict a flat SF at 
low latitude on angular scales 59 > 3' as the number of 
turbulent cells in the line of sight increases. These models 
explore angular scales much smaller than those probed 
here. We see a general increase in SF slope towards the 
Galactic plane, co ntrary to the expec tation of a smaller 
slope expected by ISun fc Reichl ()2009f ). As noted in Fig- 
ure [HI the lack of a strong longitude dependence of SF 
slope, and its wider distribution with Galactic latitude 
suggest that the slope of SFs on larger angular scales 
is affected by the local interstellar medium. This sug- 
ge sts that local struct ure not represented in the models 
of ISun fc Reichl ((2009) may be responsible for the steep 
slope of the SFs at low latitude. 

4.3. Depolarization of diffuse emission 

RMs s of extragalactic sources are affected by Faraday 
rotation along the entire line of sight from the source to 
the observer, but mostly in the Galaxy. Depolarization 
of diffuse emission extends to b « 30°. At this latitude, 
the line of sight through the Galactic ISM is approx- 
imately twice as long as the line of sight towards the 
Galactic poles. Adopting a scale height of the warm ion- 
ized medium H = 1.83±g;^| kpc (jGaensler et al.l 12008ft . 
most of the Faraday rotation in a line of sight at b = 30° 
should arise within a distance ~ 3.6 kpc. This is still 
more than an order of magnitude larger than the distance 
to the North Polar Spur, for which we adopt the dis- 
tance to the Scorpio-Cent aurus OB association (170pc; 
lEgger fc Aschenbachlll995l ). although the near side may 
be closer. Assuming an exponential density profile, and 
5n e ~ n e , SB ~ B independent of density, approximately 



12 



Stil et. al. 



10" 



Nj 100 
D 




0.1 



10 



69 (degr.) 



10- 



l 



% 100 L 



n 1 — i — i i i i 1 1 1 1 — i — i i i i 1 1 



(b) 




1 • ■ 



I i i 1 



_i i 



0.1 



10 



69 (degr.) 



Fig. 14.— SFs at the NGP (a) an d SGP (b) re g ions com- 
■aring resu l ts obt ained with data from PTavlor ct al. (200(|) and 
Gray points are SFs made with RMs from 



ao et al. 



Mao et al. I (2010), using 354 accepted RMs in the north and 319 
accepted RMs in the south. These SFs differ from Figure 6 in 
IMao et al. I 1120101 ) as discussed in the text. The bla ck dots rep- 
resent SFs that include all RMs in IMao et al~l 112010 1. The black 
curves represent the SFs from this paper (Figure [6} with the noise 
power multiplied by 1.5 before subtraction. The black curve in 
panel (a) does not extend to angular scales smaller than 0?7 be- 
cause the noise power times 1.5 is comparable to the total variance 
of the data. 



half of the Faraday rotation and RM variance occurs in 
front of the North Polar Spur. Large-scale inhomogene- 
ity of the ISM introduces significant uncertainty in this 
estimate. 

A foreground screen with RM variance erj| M on scales 
much smaller than the 0? beam size of the DRAO 26-m 
telescope depolarizes emission observed at wavelength A 
by a factor 

' "2or M A 4 ), 



— = exp(- 
Po 



(G) 



(jSokoloff et al.lll998f ). The sudden drop of a factor - 3 
in polarized intensity at |6| = 30° observed at A = 21 cm 
requires (J^ M = 280± 180 rad 2 m~ 4 , where the error rep- 
resents a factor 2 uncertainty in the depolarization factor. 
The observed total RM variance is 600 ± 100 rad 2 m~ 4 . 
Around the depolarization edge, the SF slope is between 
0.3 and 0.4 (Figure fTUj) . Assuming a typical slope 0.35 
around \b\ = 30°, we find the observed RM variance 
on the scale of the beam of the DRAO 26-m telescope 
(470 ± 100) rad 2 m -4 . Subtracting the variance re- 



lated to measurement errors er 2 



270 ± 18 rad 2 uT 



49 ± 20 rad 2 m~ 4 (Equation [3]), the RM variance at- 
tributed to the Galactic ISM along the line of sight is 
150 ± 110 rad 2 m~ 4 . Approximately half of this vari- 
ance is expected to arise behind the North Polar Spur, 
and would not depolarize emission from this region. 

The observed variance is a factor ~ 3 lower than what 
is required to depolarize the diffuse emission, but in view 
of the estimated errors, the difference is not very signif- 
icant. Within the errors, it is possibile that differential 
Faraday rotation along the line of sight contributes to 
depolarization of the diffuse emission. 

5. SUMMARY AND CONCLUSIONS 

This paper presents an analysis of structure in the 
magneto-ionized interstellar m edium probed by 37,543 
RMs derived from the NVSS ([Taylor et all [20091 ) . RM 
SFs are presented as a function of direction for 80% of 
the sky. We find that: 

1. Large-scale structures visible in RM amplitude and 
sign appear related to structures in the Local ISM, such 
as radio Loops I and II, the Gum nebula and some other 
Hll regions. We associate region C' (Table [lj with a soft 
X-ray shadow also visible in H I at Vlsr ~ — 2 km s — 1 . 

2. SFs near the Galactic poles indicate the presence 
of a Faraday screen at both Galactic poles. The vari- 
ance in RM from foregound ISM is ct 2 q M = 202 rad 2 nT 



and the NGP, and a 



ISM 



'ISM 

=419 rad 2 



m~ 4 at the SGP 



and the intrinsic variance of background sources u 2 



(120 rad 2 m -4 and 307 rad 2 m~ 4 respectively when us- 
ing Figure [TJ| . This difference is attributed to a local 
asymmetry of electron density and/or properties of the 
magnetic field with respect to the Galactic plane. 

3. The distribution of RM variance over the sky shows 
a systematic trend with Galactic longitude such that 
lines of sight through the inner Galaxy have higher RM 
variance than lines of sight through the outer Galaxy. 
This suggests that most or all of the line of sight through 
the Galaxy contributes to RM variance. The slope of the 
SF does not show the same trends, suggesting its origin 
is related to local structures that dominate RM variance 
on larger angular scales. 

4. Radio Loop III and the Orion-Eridanus bubble ap- 
pear to enhance RM variance between latitudes 30° and 
50°, and a higher SF slope is found in these areas. Several 
other curved filamentary structures are found in direction 
of the Galactic anti-center (Figure [4j . 

5. The steepest SFs (a2 ~ 1.3) are found in the direc- 
tion of the Gum nebula and the £ Oph H II regions, but 
the steep slope may be a result of a localized enhance- 
ment of the mean RM by these two H II regions. 

6. RM variance traced by compact extragalactic radio 
sources can account for depolarization of diffuse Galactic 
emission given the uncertainties, but the over-all vari- 
ance is on the low side allowed by the errors. Differential 
Faraday rotation along the line of sight may increase de- 
polarization as well. 
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